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Abstract 

^^ I The quantum H-^ integrable system is a 3D system with rational potential related to the non- 

(N 

crystallographic root system H^. It is shown that the gauge-rotated H^ Hamiltonian as well as one 

^ ; 

of the integrals, when written in terms of the invariants of the Coxeter group H-^, is in algebraic 
form: it has polynomial coefficients in front of derivatives. The Hamiltonian has infinitely-many 
O ' finite-dimensional invariant subspaces in polynomials, they form the infinite flag with the charac- 

i-G ■ teristic vector a = (1,2,3). One among possible integrals is found (of the second order) as well 

as its algebraic form. A hidden algebra of the H^ Hamiltonian is determined. It is an infinite- 
dimensional, finitely-generated algebra of differential operators possessing finite-dimensional repre- 
^ ' sentations characterized by a generalized Gauss decomposition property. A quasi-exactly-solvable 

C^ ' integrable generalization of the model is obtained. A discrete integrable model on the uniform 

o 



o 
o 



X 



lattice in a space of //3-invariants "polynomially"-isospectral to the quantum H^ model is defined. 
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I. INTRODUCTION 



About 30 years ago, Olshanetsky and Perelomov developed the Hamiltonian Reduction 



Method, later known as the Projection Method (for a review, see [jj]). This method provides 
an opportunity to construct on a regular basis the non-trivial multidimensional quantum 
(and classical) Hamiltonians, which are associated to the crystallographic root spaces of 
the classical {A^.B^.Cn.Dn) and exceptional ((72, -F4, _E6,7,8) Lie algebras. All these sys- 
tems are symmetric with respect to the corresponding Weyl group transformations. The 
Olshanetsky-Perelomov Hamiltonians have the property of complete integrability (the num- 
ber of integrals of motion in involution is equal to the dimension of the configuration space). 
There are three types of the Hamiltonians with rational, trigonometric and elliptic poten- 
tials, respectively. The Hamiltonians with rational and trigonometric potentials are exactly 
solvable (the spectrum can be found explicitly, in a form of a first- or second-degree poly- 
nomial in the quantum numbers, respectively). 

A Hamiltonian with rational potential associated to a Lie algebra g of rank A^ with root 
space A has a form 

1 ^ r r)2 -| 1 1^12 

1 ^^ — V U 2 2 \ — *" \^\ /i\ 

^^ = o2^ -T^ + ^'^k +o2^^H7— 3^' (1) 



2 

k=l 



dxl 



2 ^ "'"'(a-x)^ 



where i?+ G A is the set of positive roots, u G M"*" is a real parameter, g\a\ = i^|a|(i^|a| — 1) 
are coupling constants depending only on the root length, and x = (xi,X2, . . . ,xn) is the 
coordinate vector. The configuration space is the principal Weyl chamber of the root space 
(see ll). The ground state eigenfunction and its eigenvalue are given by 

QG-R+ 

and its eigenvalue 

^0= (f+ E ^h)^' (3) 

where ^2 is the invariant of the degree two (for definition see below). It is indicated 
in [l| that the Hamiltonian ([T]) with the property (2)- (3) can be introduced for the non- 
crystallographic root systems iJa, -^4 and dihedral l2{m) with a Coxeter group as a symmetry 
of the system. All that is not true for the H^, H4 and l2{m) Hamiltonians with trigonometric 
or elliptic potential. The complete integrability of (1) for the case of non-crystallographic 
root systems has been proven in 2| using the formalism of quantum Lax pairs. 



Following [3|, we make three definitions. 

Definition 1. A multivariate linear differential operator is said to be in algebraic form 
if its coefficients are polynomials in the independent variable(s). It is called algebraic if by 
an appropriate change of the independent variable(s), it can be written in an algebraic form. 

Definition 2. Consider a finite-dimensional (linear) space of multivariate polynomials 
defined as a linear space spanned in the following way: 

where the a's are positive integers and n G N. It represents the Newton polytope in a 
form of a rectangular pyramid. Its characteristic vector is the d-dimensional vector with 
components q;j: 

a = (ai,a2,---ad) • (4) 

For some characteristic vectors, the corresponding polynomial spaces may have a Lie- 
algebraic interpretation, in that they are the finite-dimensional representation spaces for 
some Lie algebra of differential operators. The smallest characteristic vector is Uq = 
(1, !,...,!). It corresponds to the finite-dimensional representation space for the Lie al- 
gebra gl{d + 1) of the first order differential operators acting in R*^ j^. We will call such a 



n , 



space the basic space as well as the associated flag will be called the basic flag. 

Definition 3. Take the infinite set of spaces of multivariate polynomials P„ 
n eN, defined as above, and order them by inclusion: 

Po C Pi C P2 C . . . C P„ C . . . . 

Such an object is called an infinite flag (or filtration), and is denoted P("\ If a linear 
differential operator preserves such an infinite flag, it is said to be exactly-solvable. It is 
evident that every such operator is algebraic (see j^). If the spaces P„ can be viewed as the 
finite-dimensional representation spaces of some (Lie) algebra g, then g is called the hidden 
algebra of the exactly-solvable operator. 

If a linear operator preserves several flags and among them there is a flag for which dimP„ 
is maximal for any given n, such a flag is called minimal. Every flag can be characterized 
by a normal vector n^ to the hyperplane aipi + a2P2 + . . . + adPd = n of the base of the 



Newton polytope. This normal vector is, in fact, the characteristic vector. It is clear that for 
minimal characteristic vector the angle with basic characteristic vector [a, do) is minimal. 

For any root system A there exist A^ = rank(A) homogeneous, algebraically independent 
polynomials which are invariant with respect to the Coxeter group. They are called the 
invariants. The lowest possible degrees a of these invariants are the degrees of the Cox- 
eter group. Each invariant is defined ambiguously, up to a non-linear combination of the 
invariants of the lower degrees. 

One of the ways to find an invariant of degree a (denoted as ti ) is to make averaging 
over an orbit Q, 

(see e.g. |5|), where a:'s are some formal variables which can be identified with the Cartesian 
coordinates. It is worth mentioning that for any Coxeter group there exists a second degree 
invariant ^2 ? this invariant does not depend on the chosen orbit. Later on we will use the 
invariants (Q as new variables. We will call them the orbit variables. 

For all the crystallographic root systems algebraic representations of all quantum Hamil- 
tonians, both rational and trigonometric, have been found ([6|]-[l^). The general strategy 
which was used to find a minimal flag for the rational Hamiltonians is the following: (i) as a 
first step we consider the similarity-transformed version of ([T]), namely h oc \I'g'^('H — £'0)^0, 
(ii) then, we choose a certain orbit to construct a particular set of variables which lead to 
an algebraic form of the transformed Hamiltonian h, (ii) finally, exploiting the ambiguity 
in the definition of polynomial invariants of the fixed degrees we search for variables for 
which the flag of invariant subspaces of (1) is minimal flag. A primary goal of this paper is 
to show that the same strategy can be applied for a study of the rational H3 Hamiltonian 
(1). We find the algebraic form of the rational -^3 Hamiltonian and a minimal flag of its 
invariant subspaces. Furthermore, we demonstrate that any (invariant) subspace from the 
minimal flag is a representation space of a finite-dimensional representation of a certain 
infinite-dimensional, finitely-generated algebra. This algebra is the hidden algebra of the H3 
system. A similar analysis is done for one of the integrals of the H3 system. 

Another goal of the paper is to find a quasi-exactly-solvable generalization of the H3 
rational model. By definition a linear differential operator is quasi-exactly-solvable (QES) if 
it preserves a finite-dimensional functional space with an explicitly indicated basis (see e.g. 



13|). Thus, it implies that the QES operator has a finite-dimensional invariant subspace 
spanned by known functions. Furthermore, one can indicate explicitly a basis where the 
operator being written in the matrix form has a block-triangular form. In practice, for 
all known examples of the QES operators the finite-dimensional invariant subspace is a 
space of inhomogeneous polynomials in one or several variables. In many cases the space of 
polynomials can be identified with a finite dimensional representation space of a Lie algebra 
of differential operators of the first order. In the case of crystallographic root systems a 
certain QES generalization has been found for each particular rational Hamiltonian |l4l |. 
All those examples are related with the existence of the hidden sl{2) algebra. We show that 
a similar s/(2)-quasi-exactly-solvable generalization of the if 3 model exists. 

Finally, we show that in the H3 orbit space (in the space of the if 3 invariants ti ) 
there exists a discrete model defined on three-dimensional uniform lattice with polynomial 
eigenfunctions which is isospectral to the rational H^ model and integrable. We will call 
this model the " discrete ii3 rational model" . 

II. THE HAMILTONIAN 

The Hamiltonian of the rational if 3 model (see ([1])) is invariant wrt the if 3 Coxeter 
group, which is the full symmetry group of the icosahedron. This discrete group is subgroup 



of 0(3) and its dimension is 120 (see e.g. 15|). In the Cartesian coordinates xi,X2,X3 the 
Hamiltonian has the form 

1 ^ " 



fc=l 



9^ , 22,9 



+ 00 Xl. + 



-— 2 ^u. Xfc^-2 



2a ^'^ 

where {i,j,k} = {1,2,3} and its even permutations. Here g = viy — 1) > —1/4 is the 
coupling constant, Lp± = (1 ±-\/5)/2 the golden ratio and its algebraic conjugate. We choose 
as the configuration space the fundamental domain of the if3 group - the space bounded by 
three planes 

Xi = , X3 = and X3 -|- </3+a;i -|- (p-X2 = (7) 

for xi, X2,X2, > 0. 



The ground state eigenfunction and its eigenvalue are 

3 



vl>o(x) = A^-A^ exp [ - I ^ xl) , Eo = ^u;(l + 10^) 
^ fc=i ^ 



) 



(8) 



where 

3 

Ai = Y[xk, (9) 



fc=i 



A2 = n II [xi + i-'^rv+^j + i-'^r'¥^-xk]. (10) 

{i,j,k} ^i's=o,l 
The ground state eigenfunction (^ does not vanish in the configuration space ([7]). 

The main object of our study is the gauge-rotated Haniiltonian i^ with the ground state 
eigenfunction ([8]) taken as a factor, 

hu, = -2{^or\'Hu.,-Eo){^o) , (11) 

where Eq is given by (|H]). The gauge rotated operator f ITT]) is the second-order differential 
operator without free term. By construction its lowest eigenfunction is a constant and the 
lowest eigenvalue is equal to zero. Now let us introduce new variables in ( ITT]) . 

The H2, root space is characterized by three fundamental weights Wc, c = 1, 2, 3 (see e.g. 



16|). Taking action of all group elements on fundamental weight Uc we generate orbit Qc of 



a certain length (length = T^elements of the orbit). The results are summarized as 



weight 




orbit length 


wi = (0, (/?+, 


1) 


12 


W2 = (1, <^+, 


0) 


20 


W3 = (0, 2ip+ 


0) 


30 



In order to find i^s-invariants ([5]) we choose for simplicity the shortest orbit Q{wi) and make 
averaging, 

w£ Q(uii) 

where a = 2, 6, 10 are the degrees of the H3 group. These invariants are defined ambiguously, 

up to a non-linear combination of the invariants of the lower degrees 

An) An) 

''2 '"'■ ''2 5 

if ^tf + A(tf)3, (13) 



where A, B, C are parameters. Now we can make a change of variables in the gauge-rotated 
Hamiltonian fITT]) : 

The first observation is that the transformed Hamiltonian h}i.^{t) flTTl) takes on an algebraic 
form for any value of the parameters A, B, C in variables t's f fT3|) . The second observation 
is that for any value of the parameters A, B, C there exists a flag of invariant subspaces in 
polynomials of the Hamiltonian h}i.j{t). Our goal is to find the parameters for which hii^{t) 
preserves the minimal flag. After some analysis we found such a set of parameters 

A = -1^, 5 = -^, C = i^. (14) 

10 15 ' 375 ^ ^ 

The t-variables (TT3l) for such values of parameters, which we denote as r- variables, are 

_ 2 I 2 I 2 

r2= - — {xl + xl + xl) + Y^(2 - 5(/.+) {xlxt + xjx^ + xlx^) 

+ ^^(2 - 5^_) (X1X3 + X2X1 + X3X2) - — {X1X2X3) , 

^3 = Y^ (a;}° + 4° + xf) + §(1 + 5^-) (x?x2 + xlxj + xlxl) (15) 

+ —(1 + 5^+) {xfxl + xlxl + xlxl) + —(1 - 5v9_) (x?x^ + xlxl + x^x^) 
4 112 

+ '^7Z\^ ~^ 5(/9+j (^X]^X3 + X2X1 + X3X2J rr- (,2;i2^2''^3 "r" X2X^X^ + X3X]^X2J 

, -^J-^ /244, 244, 24 4\ 

1^ ( rp rf rf _\__ rf rr> rr> _\_ rv* rf* rp \ 

"T r)f- l"''l"''2-^3 ~r "t'2-^3-^1 ~r •^3-^l"''2/ • 

The Hamiltonian hY{^{T) has infinitely-many finite-dimensional invariant subspaces 

p(i,2,3) _ (^Pv|Vf« I < pi + 2p2 + 3p3 < n) , n = 0, 1, 2, . . . . (16) 

which form the minimal flag. Its characteristic vector is 

amin = (1,2,3) . (17) 

It is worth noting that each particular space Vn flT6|) as well as the whole flag are invariant 
with respect to a weighted projective transformation 

Ti -)■ Ti + a , 

T2 ^ ^2 + 61 T^ + &2 n + ^3 , (18) 

7-3-^^3 + ci T1T2 + C2 rf + C3 r2 + C4 Ti + C5 Ti + ce . 



where {a,b,c} are parameters. It manifests a hidden invariance of the Hamiltonian ([6]). It 
is seen in a clear way in the space of orbits only. 

Finally, the gauge-rotated Hamiltonian f llip in the r-coordinates is written as 



«,J = 1 -" 4 = 1 



with the coefficient functions 



^1 = 4ri, 

Ax2 = 12r2, 

As = 20r3, 

48 2 45 

^2 = -Tir2 + —T3, 

16 o o 

15 

64 128 o 

As = -yriTsra + — r| , (20) 

Bi = 6(1 + lOi/) - 4u;ri , 

^2 = - — (l + 5i^)ri2-12a;r2, 
5 

64 
Bs = - — {2 + 5iy)nT2-20uT3. 
15 

It can be easily checked that the operator flT9|) is triangular with respect to action on 

monomials rfVl^r^^. One can find the spectrum of (1191) /iHg'/' = —2eif explicitly 

eni,n2,n3 = 2u}{ni + 3^2 + 5n3) , (21) 

where rij = 0, 1,2, . . .. Degeneracy of the spectrum is related to the number of solutions 
of the equation ni + 3n2 + bn^ = k for k = 0, 1, 2 . . . in non-negative numbers ^1,2,3- The 
spectrum e does not depend on the coupling constant g and it is equidistant. It coincides 
to the spectrum of 3D anisotropic harmonic oscillator with frequencies {2uj,6uj, IOcj). The 
energies of the original rational H-^ Hamiltonian ([6]) are E = Eq + e. 

The boundary of the configuration space of the rational H-^ model (^ in the r variables is 
determined by the zeros of the ground state eigenfunction, hence, by pre-exponential factor 
in ([8]). It is the algebraic surface of degree 15 in Cartesian coordinates being a product of 
monomials. In r-coordinates it can be written as 

12960ri'r|-5760riV2V3+640riV2^+54000riV2r|-21600rir2V3+2304r2^+50625r| = 0, (22) 



which is the algebraic surface of degree seven; the equation contains monomials of the degrees 
7, 5 and 3. It is worth mentioning that l.h.s. of fl22|) is proportional to the square of Jacobian, 



III. INTEGRAL 

The Hamiltonian ([6]) being written in spherical coordinates (r, 6, (f)) takes a very simple 
form 



n 



H-i 



_1a(3) + 1^2^2 ^ 

2 2 



where A*^^^ is the 3-D Laplacian and the angular function 

2u(u-l) 2u(u 



W{9, 



+ 
+ 



{sgC^ + (p+SgS^ + 

2u{u-l) 


V-Cgy 


2z/(z/-l) 


f-cey 


2u{u-l) 


-SgC^y 


2iy{iy- 


-ip. 

1) 


-SgC^y 


{Cg + ip+SgC^, + ip. 

2u{u-l) 


-SgS^y 



+ 



+ 



{sgC^ 


- ip+SgS^ + (p-Cgy 

2u{u-l) 


{sgC^ 


- (fi+SgS^ - (f-Cgy 

2u{v-l) 


{sgS^ 


- if+Cg + ip_SgC^y 

2z/(z/-l) 


{sgS^ 


- (f+Cg - (f-SgC^y 

2u{u-l) 


(cg- 


ip+SgC^ + Lf^SgS^y 

2z/(z/-l) 



Cg + (f+SgC^ - (f-SgS^y 

u{u — 1) u{u — 1) u{u 



Cg - (p+SgC^ - (f-SgS^y 
-1) 



2sic 



2^2 



2sis 



2^2 



2ci 



(23) 



(24) 



Here, for the sake of simplicity we denoted c^ = cos ^, s,j = sin ?9. It is seen immediately, 
that the Schroedinger equation ( l23l) admits a separation of radial variable r: any solution 
can be written in factorized form 



^{r,e,<P) = R{r)Q{e,, 



(25) 



Functions R and Q are the solutions of the equations 



1 d 



2f2 Qj, \ Qj, 



d 



2 2 

-UJ r 


+ 


7" 


2 




hJ 



R{r) = ER{r) 



J-Q(e,0)=7Q(^,0) 



(26) 

(27) 



respectively, while 7 is the constant of separation. The operator T has the form 

T = ^-C' + W{e,<jy), (28) 

where C is the angular momentum operator: 

.2 1 d f d\ 1 d^ 



sin 9 09 \ 89 J sin^ 9 d(j)'^ ' 
It can be immediately checked that the Hamiltonian T-Lhs and J-" commute, 

[Hh^,J'] = 0. (29) 

Hence, J-" is an integral of motion. Thus, it has common eigenf unctions with the Hamiltonian 

T-(-H3- 

Let us make a gauge rotation of the operator J-" (l28l) with the ground state function \l/o 
as a gauge factor, 

/ = (*o)-'(-^ - 7o)^o , 70 = Yi^(l + 15^^) , (30) 

where 70 is the lowest eigenvalue of J-", and make a change of variables to the r variables 
( IT5I) . The operator / has an algebraic form, 

ij=l J j=l J 

where 

i^ii = 0, 
F12 = 0, 

Fi3 = 0, 

24 . 45 ,„ 2 

5 4 

F23 = "4^1 ^1 + 12rfr3 + SOrars , (32) 

15 

64 3 32 2 rn 2 

F33 = - — TiTl + —r^T2T3 + 50Tl, 

Gi = 0, 

24 
G2 = —{l + 5u)Tf + 3{7 + 30u)T2, 
5 

32 

G3 = — (2 + 5z/)riV2 + 5(ll + 30z/)r3. 

10 



It is worth noting that in the operator / the variable ri appears as a parameter. It implies 
that any eigenfunction of hu^ which depends on ri only (see below Ch.V for a discussion) 
is the eigenfunction of the integral / with zero eigenvalue. 

It can be also shown that the operator / has infinitely many finite-dimensional invariant 
subspaces in polynomials 

pa,3,5) ^ (^f Vf rf^ I < j9i + 3p2 + 5p3 < n) , n = 0, 1, 2, . . . . (33) 

which form a flag with characteristic vector (1, 3, 5). 

The spectrum of the integral J-" can be found in a closed form, 

-ioMM = 2(3A;2 + hhf - 'i^h^h^ + (1 + 30i^)(3A;2 + 5^:3) + 70 , (34) 

where ki, ^3 = 0, 1, 2, . . . and 70 is given by (29). 

It can be shown that the Hamiltonian hu„ has a certain degeneracy - it preserves two 
different flags: one with minimal characteristic vector (1,2,3) and another one with charac- 
teristic vector (1,3,5). The fact that the operator hu-^ with coefficients ( !20|) commutes with 
/ given by ( l3Tll implies that common eigenfunctions of the operators hu-^ and / are elements 
of the flag of spaces p(i'3'5)_ 

Let us denote 0„.j the common eigenfunctions of hu-^ and / which are elements of the 
invariant space Pn ' and their respectful eigenvalues e„_j, 7„_j. The index i numerates these 
eigenfunction for given n starting from 0. The function </>„ j is related to the eigenfunction of 
the Hamiltonian 'Hh3 (El) (and the integral J-") through \E'„^j = \&o0n,i- Thus, the eigenfunc- 
tions {(f)} are orthogonal with the weight factor |\E'oP- As an illustration let us give explicit 
expressions for several eigenfunctions 0„ « and their respectful eigenvalues. 



n = 



n = 1 



n = 2 



^0,0 



1 , eo,o = , 70,0 = . 



3 

>i,o = ri- -—(1 + lOi^) , ei^o = 2w , 71,0 = 



5 ,. „ . 15 



^2,o = rt--il + Qi^)ri + -—{l + Qu){l + lOu), 63,0 = 4w , 72,0 = 

UJ 4:(jJ 



11 



• n = 3 

3 15 15 

03.0 =rl - —(7 + 30z/)r2 + -—{I + 6z/)(7 + 30z/)ri - --(1 + 6z/)(7 + 30i^)(l + lOz/) 

£3,0 = 6a; , 73,0 = . 

8(1 + 5z/) o 

03.1 =^2 + ,;- , ^„ ' ri , £3,1 = 6w , 73,1 = 21 + 90z/ . 

5(7 + 30z/) 

As stated before, the Hamiltonian /i^/g preserves two flags witli cliaracteristic vectors 
1^1, 2, 3) and (1, 3, 5), respectively. The angle between the normal vectors of the minimal flag 
1^1, 2, 3) and of the basic one (1, 1, 1) is given 

6 
cosOh = -^= or Oh ^ 0.39 , 
V42 



while between the vectors (1, 3, 5) and (1,1,1) 



9 
7105 



cos 6'/ = —j= or Of ^ 0.50 . 



It seems evident that if one or more extra integrals exist they will take an algebraic form 
in r-variables after the gauge rotation. 

IV. DISCRETE UNIFORM H^ SYSTEM 

The existence of the algebraic form of the H^ Hamiltonian in the space of invariants 
allows us to construct a discrete system with a remarkable property of isospectrality. This 
construction is based on employment of a quantum canonical transformation as a basis to 



perform a discretization of a continuous system [ITJ. Such a procedure was called a Lie- 
algebraic discretization. It was already used in the past to construct the isospectral discrete 
model of the harmonic oscillator (the Ai system in the Hamiltonian Reduction nomenclature) 

n 

in the space of r = x^ [l8l ]. 

Let us introduce a set of the finite-difference operators 

^i ^t (35) 

A-f ') /(r,) ^ r,/(r, - 5,) = {ne~'^^^)f{n) , 

where Si, i = 1,2,3 are spacings; here no summation over repeated indexes is implied. 
The operator V^ is the finite-difference derivative or discrete momentum; sometimes, it is 

12 



called the Norlund derivative. The operator X,^ is a discrete analogue of the multiplication 
operator. The operators V^ and XJ^ form a canonical pair, 

[I5f ),pf )] = , [Xl''\x^'^^] = , [Vf'\x^'^^] = S,, , (36) 

for i,j = 1,2,3. Hence, the operators ( l35l) span the 7-dimensional Heisenberg algebra 
realizing a three-parametric quantum canonical transformation with parameters 5i,2,3- In 
the limit when all Si tend to zero the operators (135!) gives rise to a standard coordinate- 
momentum representation. 

Take a linear differential operator C{di,Ti). Consider the eigenvalue problem 

£(a„r,)(p(r) = A<^(r), (37) 

and assume it has polynomial eigenfunctions. Performing the canonical transformation ( 135|1 
we arrive at 

£(Pf\A;(''))^(A'f))|0) = Av^(A'f))|0) (38) 

In order to make sense to this equation one should introduce the vacuum |0): 

pf'^lO) = , i = 1,2,3 . (39) 



Then the equation ( 1381) has a meaning of an operator eigenvalue problem in the Fock space 
with vacuum (l39l) . Now let us show that the eigenvalue problem (l38l) has polynomial eigen- 
functions and their eigenvalues are the same eigenvalues as for the polynomial eigenfunctions 
as the original (continuous) problem ( 137|) . 

In order to exploit the representation (l35l) let us first define the vacuum |0). The condition 
(l39l) in explicit form is 

fin + (5i, T2, Tg) = /(Ti, Ta, Tg) , 
fin, T2 + ^2, rg) = /(Ti, Ta, Tg) , 
fin, T2, T-s + ^s) = fin, T2, Ts) . 

Any periodic function with periods 6i in the coordinates Tj is the solution of these equations; 
however, without loss of generality, we can make the choice 

fin,T2,T3) = l. (40) 

13 



Let us now define the quasi-monomial 

r("+i) = T{r - (5)(r - 2(5) ■ ■ ■ (r - n6) 
Taking into account the relation 

and choosing of vacuum (140 p , it is easy to check that 



(41) 



X, 



{S^) 



y\o)=rt' 



(42) 



Now we can relate the solutions of (138|) with the solutions of ( 137|) . Let us assume that 



(43) 



is a polynomial solution of the equation fl37|) . The canonical transformation fl35l) implies 
the replacement of Tj by X^ . Taking in account (l42l) we come to the conclusion that 
each monomial in f H3|) should be replaced by a quasi-monomial. Hence, the corresponding 
polynomial solution of ( 155|1 is 

^(r) = 5^ aum rf Vi'Vf ) , (44) 

with the same expansion coefficients akim as in f H3|) and the same eigenvalue. 

Performing the procedure of canonical discretization (137|) — ;■ (!38|) for the /fs Hamilto- 
nian in the algebraic form hn-j fITU]) . we arrive at the following (isospectral) finite-difference 
operator: 






(45) 



with the following non-vanishing coefficients 



^0,0,0 
^1,0,0 
^-1,0,0 
^-2,0,0 
^0,-1,0 
^0,0,-1 



-^(2 + 5ia;) 

0\ 



2 

Ji 
4 

^\ 
48 

5^ 
12 

(51^2 

5 



n 3T2 573 

(5i (^2 ^3 

2ri 12r2 20r3 



-^(1 + 10^) 



5o 



5, 



3(1 + lOi/) 



^(l + 5ia;)ri 



t-i(t-i - h 



2t2 , 5r3 



^2 



^3 



l + hu 



(2 + 5ia;)r2 , 



.(2 + M + ^ 



5l5; 



-^3 



(46) 
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^0,-3,0 

^1,-1,0 

^1,0,-1 

^-1,-1,0 

^-1,-2,0 

^-1,-1,-1 

^-1,-2,1 

^-2,1,0 

^-2,-1,0 

^-2,0,-1 

^-2,1,-1 

^,1,-1 

^0,2,-1 

^0,-3,1 

^0,-3,2 



128 
45^ 



T"2(t-2 -52)(t-2 -2^2) 



2472 

^1^2 
40X3 



6163 

32 

"155; 
32 

1552^3 

32 



riT2 



't2 20r3 



S2 S3 

TiT2{t2 -62) , 
T2 10X3 



m + 2u) 



riT2 



S2 



15^3 

64 



S3 



5(1 + 2z/) 



^3 
32 



1552^3 



■nT-2(r2 -^2) 



-55-2"'^^'-^'' 

772'^2Ti(ri -^1) 



T2 573 

<^2 ^3 



l + 5z/ 



5S2 

48 



^2^3 

48 



r3T-i (n - 5i 



(^2^3 

_4573 

r2 ' 
"^2 

4573 

2(^2 ' 

256 



Tsnin -61) 



4561 

128 
45^ 



7-2(^2 -S2){t2 -262) , 

7"2(r2 -^2)(t-2 -2^2) • 



The corresponding eigenvalue problem is 

y^ ^fcifefcs '^i'Tl + ^1^1' ^2 + /^2<^2, ^3 + ^3*^3 
ki,k2,k3 

which is in the explicit form 



-2e ip{Ti,T2,T3) 



2 

+ 



(2 + Situ) 



Ti , 3r2 , 5r3 



61 



61 62 S3 
2ti 12t2 2OT3 



3(l + 10i^) 



^iTi,T2,T3) 



5o 



63 



3(l + 10z/) 



if{Ti +6i,T2,T3) 



+ ^(1 + 5iUj)ti v9(ri - (5i, r2, r3) + ^^(2 + 5iUj)t2 ^{ti, T2 - 62, T3) 



12 
6i52 



(47) 



(48) 



15 



+ 3-ri(ri-5i) 

5(32 



272 ^573 

5o 5?. 



:i + 5z/) 



V9(ri -25i,r2,r3) 



+ 



2 

128 



5l(5; 



■(2 + (5iu;) + -2 

On 



^3 </'(n,r2,r3 -(^s) 



2472 



ip{Ti + (5i,r2 -(52,r3) 



,^.2"'^2(i"2 - '52)(t-2 - 2(^2) (/?(ri,r2 - 3^2,^-3) 
45()3 did; 



V9(ri + di,r2,r3-d3) 



+ 



32 

32 



riT2 



T2 20r3 



5(1 + 2z/) 



15(525. 



T\T2{r2 - 82) ^{ri - 5i,T2- 262, T3) 



C/?(ri -6i,T2 -S2,T3) 



2^3 



32 

15^ 
64 



-7-1X2 



7-2 10r3 



^2 ^3 



5a + 2z/) 



V9(ri -(5i,r2 -52,r3 + 53) 



~ ^^1^27-3 (^(ri - 5i, T2 - 82, T3 - S3) 



'3 
32 



I562S3 



-tiT2{t2 - 62) (fin - 6i,T2- 262, n + 63) 



48 r A^ 



^2 5r3 

^ + ^ + (l + 5i^) 

^2 O3 



ip{Ti -2(5i,r2 + (52,r3) 



32 



-nT-2(T-2 - ^2) Vin - 61, T2 - 262, T3 + 63) 



:i + 5u) 



1562S3 

- 3^7-1 (ri-5i) 

5^2 

48 

- ■rT2T27-l(ri - 61) (p{ti - 26i, T2 - 62, T3] 

D02 



72 573 

So S3 



ip{Ti - 2Si,T2 + S2,T3] 



48 

^2^3 

48 

S2S3 
45r3 



T37-i(ri - Si) Lp{ri - 2Si, T2, T3 - S3) 

7-37-l(ri - Si) ip{Ti - 2Si, T2 + S2, Ts - S3) 



-p- ^in, T-2 + ^2, T3 - S3) + —2" </?(ri, T2 + 2S2, T3 - S3) 

256 



4573 

■2 
'2 



4551 
128 
455| 



7-2(7-2 - 52)(r2 - 2S2) ^iri,T2 - 3S2, T3 + S3) 

7-2(7-2 - 52)(r2 - 2^2) f{Ti,T2 - 3^2, T3 + 2^3) = -2e ip{Ti, T2, T3) . 



It defines the discrete uniform H3 system. It is worth noting that, although we started 
from a second-order differential operator, the 22-point finite-difference operator occurs - it 
connects the function in 22 different points in the lattice space: four points in ri-direction, 
six points in r2-direction, four points in r3-direction. The structure of the operator is shown 
in FigUl 

The spectrum of the discrete operator hn^ ( 1451) for polynomial eigenf unctions coincides 
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FIG. 1: Graphical representation of the 22-point discrete operator (j45p . The origin is shown in 
red. 



to the spectrum of the continuous operator hii^ (fT9|) where all eigenf unctions are polynomial 
ones: 

eni,n2,n3 = 2a;(ni + 3?22 + Sns) , (49) 

where the n's are non-negative integers. The eigenf unctions of (H5|) are related to the eigen- 
functions of the continuous operator hn^ by replacing each monomial with a quasi-monomial 
in each variable. Such a phenomenon can be called partial or polynomial isospectrality. We 
cannot exclude the existence of other eigenstates of the discrete operator hn^ than those 
given by polynomial eigenfunctions. These eigenstates can correspond to non-normalizable 
eigenfunctions of hff^ . 

As a particular case let us consider the unit spacing Si = 62 = S3 = 1. This case 
corresponds to the discretization in a cubic lattice the space of orbits (in r-space) with unit 
lattice vector. The equation ( 148|) is reduced to the equation 

2 [(2 + u) (n + 3t2 + Sra) - 3(1 + lOz/)] ^{n, T2, r^) 
+ 2 [2ri + 12t2 + 2OT3 + 3(1 + lOz/)] ^(n + 1, ra, rg) 
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+ 4(1 + uj)ti ifiin - 1, r2, Ts) + 12(2 + uj)t2 ^(n, ra - 1, Ts) 

5 
+ -(25 + 8u;)r3 (/^(ri, ra, rg - 1) - 24r2 </?(ri + 1, ra - 1, rg) 

48 
+ ^n(ri - 1) [2t2 + 5r3 + (1 + 5z/)] v?(ri - 2, ra, rg) 
5 

1 28 

T2{t2 - 1){t2 - 2) V9(ri, r2 - 3, rg) - 40r3 (f{Ti + 1, r2, Ts - 1) 



45 



32 

- Y^riTa [r2 - 20r3 - 5(1 + 2z/)] V9(ri - 1, ra - 1, rg) 

32 48 

+ 7rrir2(r2 - 1) (^(ri - 1, ra - 2, rg) - —T2Ti{ti - 1) (/)(ri - 2, r2 - 1, Tg) 
15 5 

32 

+ Y^nT2 [t2 - lOrg - 5(1 + 2v)\ (^(n - 1, r2 - 1, rs + 1) 

- ^T-ir2r3 V3(ri - 1, r2 - 1, rg - 1) - —TiT2{t2 - 1) v^(ri - 1, r2 - 2, rg + 1) 

3 15 

48 

- ^T-i(n - 1) [7-2 + 5r3 + (1 + 5z/)] (^(ri - 2, r2 + 1, rg) 

5 

- 48r3ri(ri - 1) v?(ri - 2, r2, rg - 1) + 48r3ri(ri - 1) (/^(ri - 2, r2 + 1, n - 1) 

- 45r3 (/?(ri, r2 + 1, r3 - 1) + — ^ (/^(ri, r2 + 2, r3 - 1) 

256 128 

+ ^F^2(r2 - l)(r2 - 2) ^(n, r2 - 3, r3 + 1) + ^^r2(r2 - l)(r2 - 2) (^(n, r2 - 3, rg + 2) 
45 45 

128 , ,, , , 

+ -rz-'T2[T2 - 1)(t-2 - 2) V5(ri, T2 - 3, r3 + 2) = -2e (/)(ri, r2, t^,) . 
45 



(50) 

A similar procedure of discretization can be applied to the integral J-". Instead of the 
continuous algebraic operator / we get its discrete counterpart 



{S^) v(^■>)^ 



with the following coefficients 

-So, -1,0 

Bo-2fi 

Bo,o,-i 

-Bo,0,-2 
-S-1,1,0 



fci,fc2,fc3 



fci(5i9i+fc2<5292+fc3^393 



(51) 



3r2 



10r3 12r2 



6. 



62 



5 + 30i^ 



-p2 ^2(7-2 - S2 
O2 



5t3 



6r2 ^ 20r3 ' 

-^ H — r 9 + 30z/ 

^2 03 



50 

5: 



2^-3 (7-3-53) 



^3 

1273 
^2^3 



ri(ri-(5i)(ri-25i) 



(52) 
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-6-1,2,-1 

-B-1,0,-1 

-S-1,-3,0 
-6-1-3,1 
5-1,-3,2 
-6-2,-1,0 
-6-2,-2,0 
-6-2,-1,1 
-6-2,-1,-1 
-6-2,-2,1 
-6-3,1,0 
-6-3,-1,0 
-60,-1,-1 



3Ti73 
2 62 

45 



^(ri-<5i)(ri-2<5i) + i^ 
^3 02 



WJ'^' 



46s 
64 



nT-3 



15 16 



{n-6^){n-26,) 



4561 



riT2{T2 -S2){r2 -262) 



128 
45l 



rir2(r2 -52)(r2 -2^2) 



64 



455| 



T-iT-2(r2 -52)(r2 -252) , 



156: 



-riin - 6i)t2 



T2 40X3 



62 63 
nin - 6i)T2{r2 - 62) , 

T-2 20r3 



9-20z/ 



15^2^3 

8 

-ri(ri - 6i)t2 



156. 



62 6s 



9-20z/ 



^ri(ri-5i)r2r3 



Tlin -6i)t2{t2 -62) 

n{n-6i){n-26i) J 

62 



156263 
24 

24 
^ri(ri-(5i)(ri-2(5i)r2 

30 

^2^3 



l + 5i/ 



-T-2T-3 • 



We end up with the 22-point finite-difference operator occurs - it connects the function in 
22 different points in the lattice space: four points in ri-direction, six points in r2-direction, 
five points in r3-direction. The structure of the operator is shown in Figj2j It is surprising 
that both discrete operators hn^ and / have the same structure connecting 22 points. 

It is evident that the discrete operators hn^ and / continue to commute. Such a procedure 
of discretization preserves the property of integrability. 
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-352 -262 -62 

T2 



262 



FIG. 2: Graphical representation of the 22-point discrete operator (I5ip . The origin is shown in 
red. 

V. QUASI-EXACTLY-SOLVABLE GENERALIZATION 



Among the eigenf unctions of the Hamiltonian hu.^ flT9|) there is an infinite family of 
eigenfunctions depending on the single variable ri. These eigenstates are solutions of the 
eigenvalue problem 



9V 



dip 



-hi^ = -4ri--^ + {Aion - 6(1 + 10z/))tt^ = e^p . 



(53) 



Corresponding eigenfunctions are given by the Laguerre polynomials and the eigenvalues are 
linear in quantum number 



ipnAri) = 4^/2^0;^^) ^ ^^^ ^ 4^^^ ^ ^^ ^ Q^ ^^ 2, 



(54) 
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The operator in the Lh.s. of ( 153|) can be rewritten in terms of the generators J°, J of the 
Cartan subalgebra of the algebra s/(2) of the first order differential operators: 

A--f^-^-n. A°^nA-|, r^l^, (55) 

(see [ij])- For integer k the generators fl55|) have a common invariant subspace in polynomials 
of degree not higher than /c, 

Vk = {rl\Q<p<k) , (56) 

where dimP^ = [k + 1). The operator flH5]) takes the s/(2)-Lie-algebraic form 

/ii = 4J°J~-4a;Jo° + 6(l + 10i^)J" . (57) 

It is easy to check that the operator hi preserves an infinite flag of spaces of polynomials 

iSD, 

n C Pi C P2 C ■ ■ • C n C ■ ■ ■ , (58) 

and, in particular, any eigenfunction is an element of the flag. 

Let us proceed to construct a QES generalization of ([6]). We look for the QES Hamiltonian 
in a certain form 

H^;if = HH, + V^'^''\n), (59) 

where V^'''^^'^ is a potential. Let us make a gauge rotation of (!59|) of the form (fTTj) . We impose 
the requirement that the resulting operator possesses a ri-depending family of eigenfunc- 
tions. We obtain the following equation: 

_/,(^-)y, = _4^^^ + (4,^^^_6(l+10i.))-^ + 2V(^^^)(ri)</p = e^ . (60) 

C/7~-| uT\ 

Our aim is to find V^'^^^'^ for which the operator h^'^ is s/(2)-Lie-algebraic - can be 



rewritten in terms of the generators fl55l) . Following [IJ], let us gauge rotate the operator 

iQD, 

hf^'^-^-^ = rr-exp [Irl) /.j-V^ exp [-^rf) 

= 4ri— ^ - 2{2aTl + 2u;ri - 3 - 47 - 30z/)— + aVj^ + 2aa;ri^ (61) 

- 2a (27 + 15. + g ,, + ^^P^'+^^'"- + ') - 4c.o^ - 2V<-)(rO . 
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If the corresponding potential V^^^es) jg chosen of the form 



y(qes) ^ 1 ^2 3 ^ 

2 ^ 



1 a ( 2A; + 27 + 15z/ + - ) Ti + 



7(27 + 30z/ + l) 



n 



the operator /i^""'^'' ^'''''' has the Lie-algebraic form 



h\ 



(sl{2)-qes) 
1 

{sl(2)-qes) 



4J°J- - 4aJ+ - 4a;J° + 2(A; + 47 + 3(1 + lOz/)) J" , 



(62) 



(63) 



where a > and 7 are parameters, here the constant terms are dropped off. 

It can be seen that the operator hi '^^^ (see (l63l) ) has the space Vk as the invariant 
subspace, but it does not preserve the flag of spaces fl58|) . Hence, (1571) has (k+1) polynomial 
eigenf unctions of the form of polynomials of the degree k, 

^^'^(n) = E7pVi, j = 0,l,2,... , 



i=0 



while other eigenf unctions are not polynomials. Now we can give the final expression of the 
s/(2)-quasi-exactly-solvable Hamiltonian associated with the root space H^: 

Q2 



H 



(qes) 
H3 



1 ^ 

fc=l 



dxl 



, 2 2,9 

+ ^ ^fc + ^2 



2^ 



+ 



- a2(x2)3 + aujipi^f - a ( 2A; + 27 + 15z/ + - ] x^ 
7(27 + 30i^ + 1) 



(64) 



where {i, j, /c} = {1, 2, 3} and its even permutations, and x^ = X]i=i ^1- -^^^ ^'^is Hamiltonian 
we know (/c + 1) eigenstates explicitly. Their eigenf unctions are of the form 

^k{x) = A^'A^ (x^)^ . P,,(x2) e-^^'-t(-')' (65) 

where Pk is a polynomial of degree k, the coupling constant g = v{u — 1) > —\ and Ai^2 are 
given by dH]). It is worth presenting several Pk explicitly. 



Po(x2) = 1 , EQ = -uj(l + lQu + --i 



PiA^' 



x" + — , 
2a L 



u ± ^/u^TM4^iT3{lTlO^ 



(66) 



Ei,± = Eo + coT V^^ + 2a(47 + 3(1 + lOz/)) . 
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The solutions ioT k = 1 are related through analytic continuation in one of the parameters 
u, a, 7, z/ keeping other parameters fixed. They form two-sheeted Riemann surface. 

The QES Hamiltonian f l64p is integrable - the integral J-" f l28p remains to commute with 
the Hamiltonian. 

VI. HIDDEN ALGEBRA 

We have shown that the Hamiltonian in the algebraic form flT^ acts on the finite- 
dimensional spaces of multivariate polynomials "Pn ' , n = 0, 1,2, . . . (see (ITB]) ). A goal 
of this Section is to show that each one of these subspaces is a representation space of an 
infinite-dimensional algebra of differential operators which we call /i*-^^ and to study this 
algebra. 

The algebra m^' is infinite-dimensional but finitely-generated. Their generating elements 
can be split into two classes. The first class of generators (lowering and Cartan operators) 
act in Vn for any n E N and therefore they preserve the fiag P^^'^'^). The second class 
operators (raising operators) act on the space "PA ' for a certain value of n only; they do 
not act on a space at other n's. 

Let us introduce the following notation for the derivatives: 

fj — fj — ^ fj — ^ 

OTi OTiOTj OTiOTjOTk 

The first class of generating elements consist of the 22 generators where 13 of them are the 
first order operators 



4'^ = d, , 


7.(2) 


= ^2, 


4'^ = ds , 


T['' = nd, , 


t(2) 
-'2 


= ^2^2 , 


Tt = T,d,, 


Tf ) = nds , 


t(3) 
-^11 


= r',d, , 


7^(3) _ ,3 
^ 111 - T^i C3 , 


Tf ) = nd, , 


t(2) 
-'ll 


= r',d, , 


ri'^ = r2C?3, 




-^12 


= riT2d3 


5 


6 are of the second order 








T^' = r29n , 


7^(13) 
-^22 


= rid,, , 


-^ 222 — ^2 ^^33 ! 


Tr=r35,,, 


7^(22) 
-^3 


= T3<922 ; 


Tn(22) 
-^13 = nr3022 



(67) 



(68) 



23 



and 2 are of the third order 

Tr^=r39ni, Tir^=rld,,,. (69) 

The generators of the second class consist of 8 operators where 1 of them is of the first 
order 

Jt = n Jo , (70) 

4 are of the second order 

and 3 are of the third order 

4^11 = ^3^11 Jo , J3% = r39iJo(Jo + l) , J3+ = r3Jo(Jo + l)(Jo + 2) , (72) 

where we have introduced the diagonal operator 

Jo = ndi + 2T2d2 + 3x3^3 - n . (73) 

for a convenience. In fact, this operator is identity operator, it is of the zeroth order and, 
hence, it belongs to the first class. 

Before to proceed to study the commutation relations between generators we introduce 
a notion of conjugation. Let Ti and T2 be operators acting on a monomial. We say that T2 
is a conjugate to Ti if the operator T2T1 leaves the monomial unchanged. There is a certain 
ambiguity related with the central operator Jq (!73|) . Formally, it seems self-conjugated. 
From another side, the operator Jq is, in fact, the unit operator. Thus, one can define a 
conjugate to Jq to be equal to 1 and visa versa: a conjugate to 1 is equal Jq. /,From this point 
of view any operator is defined up to a multiplicative factor of Jq. We resolve this ambiguity 
by defining the generators fl67j) - fl73j) in a way for supporting below-presented commutation 
relations and eventually a structure of algebra. In particular, Tq is conjugated to J^*", 



di < — > Ti Jo 



and Tq is conjugated to J^, 



a3^-^r3Jo(Jo + l)(Jo + 2) . 
Except for this ambiguity, the conjugation coincides with the Fourier transform. 
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A certain number of generating operators ( !67|) - (172|) span ten Abelian subalgebras: 



T _ (rpi'i) rp{3) rp(3) rp(3) T 

-^ — l-'o '-'l '-^ll '-'ill/ 

R={Tl^\Ti^\Tif} 



{T, 



(3) ^(3) 
2 5-^12 



} 



(22) rpi22)^ 
13 . 

{Tif} 



G 



^ -^~l-^3 5 ^^3,-11' "^3,-1' ^^3 J 
-)• Jt = |i2 , <^2,-l' "^2 I 

->^ = {rf),j3%} 

-7- £ = |i22 , '-'22, -3/ 



(74) 



The remaining generators span a non-commutative algebra isomorphic to gl{2) ®TiS^\ where 
TlP'"' is the Abehan algebra of dimension 2, 



.(1) ^(1) ^(2) ^(3) 



B = {T','\Tl'\T','\T^'\ J,, Jt} 



(75) 



The arrow in (171|l means that the corresponding operators are related by conjugation (for 
example, Tq and 73*"). The subalgebra B is the unique algebra those operators are self- 
conjugated. 

Decomposition (17^ . fl75l) allows us to give a compact representation of the commutation 
relations relating the generating elements. As first we can show that the commutators 
between the Abelian algebras from the l.h.s. (r.h.s.) of (1711) are closed - they are either zero 
or written in terms of generators of one of the algebras from l.h.s. (r.h.s.) plus from the 
algebra B 

[£,€]=P2(^), 

[£, 0] = 0, 
[^, e] = 0, 

[0^,6]=P2(5), 

[^,e] = P2(^©5), 

[e,0]=P3(5©i?), 



X,P 


= 0, 


[L.F] 


= 0, 


L,E 


= ^2(P), 


L,G 


= 0, 


R,F 


= L, 


[R,E] 


= 0, 


R,G] 


= P2{Fl 


'F,E 


= P2{R®B] 


P,G] 


= 0, 


E,G 


= P-i{F®B 
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Here Pk{Q) means that the commutator is given by a polynomial of the k-i\i degree in 
the generators of Q. It turns out all commutators are symmetric under conjugation. This 
property also holds for cross commutators of the algebras from the l.h.s. and the r.h.s. in 

dZlD, 

[L, % = P^iF © B), [£, R] =P^{^(BB), 

[L, 5] = P2{R © B), [£, F] = P2{m © B), 

[L,(5]=P2iR®E), [£,G] = P2(^©S), 

[R,d]=E, [9^,F] = i£, 

[R, €] = P2{F © P), [9=1, E] = P^id © B), 

[R,e] = 0, [5=l,G] = 0, 

[F,^]=G, [d,E] = (3, 

[F, 0] = P2{E © B), [5, G] = P2{(B © B), 

[E,&]=0, [e,G] = 0, 

The commutator between any Abelian algebra from ( FM|) and the algebra B is of the type 
of a semi direct product: 

[L,B] = L, [R,B]=R, [F,B]=F, [E,B] = E, [G,B] = G, 

[il,B] = £, [D\,B]=m, [d,B] = d, [^,B] = (£, [<3,B] = &, 

At last, let us indicate commutators between conjugated algebras in (TMj) - all of them are 
polynomials in generators of B: 

[L,£] = Ps{B), [R,m\ = P2{B), [F,d] = P2{B), 

[E,(B]=P3{B), [G,(3] = P,{B). 



Latter two types of relations are represented by triangular diagrams, see for example. Fig. 3. 
In general, the commutation relations between two operators are characterized by a non- 
linear combination of the generators in r.h.s.. Calculating double, triple, etc. commutators 
one can see that the commutation relations can not be closed at any order and the order of 
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B 

L- ^ 



FIG. 3: Triangular diagram relating the subalgebras L, £ and B. It is a generalization of Gauss 
decomposition for semi-simple algebras. 

the r.h.s. is increasing. Hence, the h^^"* algebra is not a polynomial algebra. It is the infinite 
dimensional algebra of ordered monomials in the 30 generating operators (l67j) - (!72l) shown 
above. 

Since h^^^ is the algebra of differential operators acting on Pn ' it should be possible to 
write the hH3 Hamiltonian (TT9|) as a combination of the (flag-preserving) generating elements 
(l67|) - (l69|) of h^^\ The /i^'^^-algebraic form of the H^ model ([6]) is the following: 

hn, = ATi'^Tl,'^ + 24rf )r(^) + 40rf ^T^^^ - -T?^T?;' + -T^^ + -T^^T^^ 

- 48Tf )t(? - ^-^Tl'h^^ + ^tS^ + (6 + mv)Ti^^ - AujT^^ (76) 

- -(1 + 5z.)t(? - 12c^Tf - ^(2 + hu)T[f - 20a;Tf . 

10 

VII. CONCLUSIONS 

We have shown that the H^ rational system related to the non-crystallographic root sys- 



tem if 3 is exactly solvable with the characteristic vector (1, 2, 3) 27|]. This work complements 



the previous studies of t 



.g. 13 - [nl). 



le rational (and trigonometric) models, related with crystallographic 



root systems (e.g. [7| - [111)- A certain signiflcance of exploration of the H^ rational system 
is due to a fact that this model is deflned in three-dimensional Euclidian (physical) space. 
There are very few known exactly-solvable systems in this space - the Coulomb problem, 
four-body Calogero-Sutherland (^3) and -BC3 rational-trigonometric models among them. 
Surprisingly, all of them are integrable while the Coulomb, ^3 and -BC3 rational problems 
are superintegrable. The same is correct for all known two-dimensional exactly- solvable 
problems in the Euclidean space: all of them are integrable. 

27 



Taking Coxeter invariants of H^ as coordinates provided us a way to reduce the rational H^ 
Hamiltonian to algebraic form. It gave us a chance to find the eigenf unctions of the rational 
H^ Hamiltonian which are proportional to polynomials in these invariant coordinates. It 
seems correct that these eigenfunctions exhaust all eigenf unctions in the Hilbert space. It is 
worth noting that the matrix Aij{T) which appears in front of the second derivatives after 
changing variables in Laplacian from Cartesian to the H^ Coxeter invariant coordinates (see 
Eqs. fl20l) ) has polynomial entries corresponding to flat space metric, hence the Riemann 
tensor vanishes. We call metric the Arnold metric 



It should be stressed that it was stated that the Hamiltonian of the H^ rational system ([T]) 
is completely integrable [2| . This implies the existence of two mutually-commuted operators 
(the 'higher HamiltoniansO which commute with the Hamiltonian forming a commutative 
algebra. It is known (see [l|) for the crystallographic systems that these higher Hamiltonians 
are the differential operators of the degrees which coincide to the minimal degrees of the 
root space (the Lie algebra) or their squares for the A^ case. It may suggest that for the 
H^ rational system the commuting integrals might be differential operators of the orders 
six and ten. Their explicit forms are not known so far. It is evident that these commuting 
operators should take on an algebraic form after a gauge rotation (with the ground state 
function as a gauge factor), and a change of variables from Cartesian coordinates to the 
Coxeter invariant variables r's. Following the experience with different integrable systems, 
it seems the integral(s) related with separation of variables do not enter to the commutative 
algebra. Therefore, the integral J-" is out of the commutative algebra of integrals. It might 
serve as an indication to a superintegrability of the H-^ rational system. 

An analysis similar to the analysis of this paper has not yet been presented for the case of 
the rational systems related to the non-crystallographic root spaces H^. A study in progress 
indicates that the characteristic vector for the quantum H^ integrable system is (1, 5, 8, 12). 
'n the case of dihedral group l2{k) the rational model has the characteristic vector (1,/c) 



2l| . It should be pointed out that unlike the rational models it is not possible to construct 
integrable (and exactly-solvable) trigonometric systems related to the non-crystallographic 
root spaces as a natural generalization of the Hamiltonian Reduction Method [l|. 

The existence of algebraic form of the H^ rational Olshanetsky-Perelomov Hamiltonian 
makes possible the study of their polynomial perturbations which are invariant wrt the H^, 
Coxeter group by purely algebraic means: one can develop a perturbation theory in which 

28 



all corrections are found by linear algebra methods [22|. In particular, it gives a chance to 
calculate the H^, Coxeter-invariant, polynomial correlation functions by algebraic means. 

Another important property of the existence of algebraic form of the H^ rational Hamil- 
tonian is a chance to perform a canonical, Lie-algebraic discretization to uniform (see Ch.IV) 



and exponential [23| lattices. In the case of both lattices such a discretization preserves a 
property of integrability, polynomiality of the eigenfunctions remains and it is isospectral. 
Although it does not give a hint how to introduce a scalar product for a discrete model. Mak- 
ing the weighted projective transformation ( TT8|) of the i^s algebraic form ( TT9|) we arrive at 
different algebraic form of the H^ Hamiltonian. Making then the Lie-algebraic discretization 
we arrive at a discrete model related to an original discrete model via change of variables. It 
can be considered as a definition of a polynomial change of variables for discrete operators. 
We found the s/(2)-quasi-exactly-solvable generalization of the if 3 model which remains 
integrable. This is the first example of quasi-exact-solvability related to non-crystallographic 
root systems. It complements the results obtained previously for all rational models related 
to crystallographic systems (see [iJ]) and for the l2{k) rational model 2l| - each of these 



models admit a certain s/(2)-QES generalization. 

Thank to the explicit knowledge of the ground state function ([8]) supersymmetric H^ 



model can be constructed following a procedure realized in 2J] for A^ rational model, in 



25| for the BCj^ rational model and in 26| for the l2{k) rational model. It can be done 



elsewhere. 
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